Mapping the Spatial Distribution of Winter Crops at Sub-Pixel Level Using AVHRR NDVI Time Series and Neural Nets

For large areas, it is difficult to assess the spatial distribution and inter-annual variation of crop acreages through field surveys. Such information, however, is of great value for governments, land managers, planning authorities, commodity traders and environmental scientists. Time series of coarse resolution imagery offer the advantage of global coverage at low costs, and are therefore suitable for large-scale crop type mapping. Due to their coarse spatial resolution, however, the problem of mixed pixels has to be addressed. Traditional hard classification approaches cannot be applied because of sub-pixel heterogeneity. We evaluate neural networks as a modeling tool for sub-pixel crop acreage estimation. The proposed methodology is based on the assumption that different cover type proportions within coarse pixels prompt changes in time profiles of remotely sensed vegetation indices like the Normalized Difference Vegetation Index (NDVI). Neural networks can learn the relation between temporal NDVI signatures and the sought crop acreage information. This learning step permits a non-linear unmixing of the temporal information provided by coarse resolution satellite sensors. For assessing the feasibility and accuracy of the approach, a study region in central Italy (Tuscany) was selected. The task consisted of mapping the spatial distribution of winter crops abundances within 1 km AVHRR pixels between 1988 and 2001. Reference crop acreage information for network training and validation was derived from high resolution Thematic Mapper/Enhanced Thematic Mapper (TM/ETM+) images and official agricultural statistics. Encouraging results were obtained demonstrating the potential of the proposed approach. For example, the spatial distribution of winter crop acreage at sub-pixel level was mapped with a cross-validated coefficient of determination of 0.8 with respect to the reference information from high resolution imagery. For the eight years for which reference information was available, the root mean squared error (RMSE) of winter crop acreage at sub-pixel level was 10%. When combined with current and future sensors, such as MODIS and Sentinel-3, the unmixing of AVHRR data can help in the building of an extended time series of crop distributions and cropping patterns dating back to the 80s.JRC.H.4-Monitoring Agricultural Resource

Estimating Sub-pixel to Regional Winter Crop Areas using Neural Nets

The current work aimed at testing a methodology which can be applied to low spatial resolution satellite data to assess inter-annual crop area variations on sub-pixel to regional scales. The methodology is based on the assumption that within mixed pixels land cover variations are reflected by changes in the related hyper-temporal profiles of the Normalised Difference Vegetation Index (NDVI). We evaluated if changes in the fractional winter crop coverage are reflected in changing shapes of annual NDVI profiles and can be detected by using neural networks. The neural nets were trained on reference data obtained from high resolution Landsat TM/ETM images and additional ancillary data readily available (CORINE land cover). The proposed methodology was applied in a study region in central Italy to estimate winter crop areas between 1988 and 2002 from 1 km resolution NOAA-AVHRR profiles. The accuracy of the estimates was assessed by comparison to official agricultural statistics using a bootstrap approach. The method showed promise for estimating crop area variation on sub-pixel (cross-validated R2 between 0.7 and 0.8) to regional scales (normalized RMSE: 10%) and proved to have a significantly higher forecast capability than other methods used previously for the same study area.JRC.DG.G.3-Monitoring agricultural resource

Retrieval of Vegetation Biochemicals Using a Radiative Transfer Model and Hyperspectral data

Accurate quantitative estimation of vegetation biochemical characteristics is necessary for a large variety of agricultural and ecological applications. The advent of hyperspectral remote sensing has offered possibilities for measuring specific vegetation variables that were difficult to measure using conventional multi-spectral sensors. In this study, the potential of biophysical modelling to predict leaf and canopy chlorophyll contents in a heterogeneous grassland is investigated. The well-known PROSAIL model was inverted with HyMap measurements by means of a look-up table (LUT). HyMap images along with simultaneous in situ measurements of chlorophyll content were acquired over a National Park. We tested the impact of using multiple solutions and spectral sub-setting on parameter retrieval. To assess the performance of the model inversion, the RMSE and R2 between independent in situ measurements and estimated parameters were used. The results of the study demonstrated that inversion of the PROSAIL model yield higher accuracies for Canopy chlorophyll content, in comparison to Leaf chlorophyll content (R2=0.84, RMSE=0.24). Further a careful selection of spectral subset, which comprised the development of a new method to subset the spectral data, proved to contain sufficient information for a successful model inversion. Consequently, it increased the estimation accuracy of investigated parameters (R2=0.87, RMSE=0.22). Our results confirm the potential of model inversion for estimating vegetation biochemical parameters using hyperspectral measurements.JRC.DG.G.3-Monitoring agricultural resource

ARCTIS — A MATLAB® toolbox for archaeological imaging spectroscopy

Imaging spectroscopy acquires imagery in hundreds or more narrow contiguous spectral bands. This offers unprecedented information for archaeological research. To extract the maximum of useful archaeological information from it, however, a number of problems have to be solved. Major problems relate to data redundancy and the visualization of the large amount of data. This makes data mining approaches necessary, as well as efficient data visualization tools. Additional problems relate to data quality. Indeed, the upwelling electromagnetic radiation is recorded in small spectral bands that are only about ten nanometers wide. The signal received by the sensor is, thus quite low compared to sensor noise and possible atmospheric perturbations. The often small, instantaneous field of view (IFOV)—essential for archaeologically relevant imaging spectrometer datasets—further limits the useful signal stemming from the ground. The combination of both effects makes radiometric smoothing techniques mandatory. The present study details the functionality of a MATLAB®-based toolbox, called ARCTIS (ARChaeological Toolbox for Imaging Spectroscopy), for filtering, enhancing, analyzing, and visualizing imaging spectrometer datasets. The toolbox addresses the above-mentioned problems. Its Graphical User Interface (GUI) is designed to allow non-experts in remote sensing to extract a wealth of information from imaging spectroscopy for archaeological research. ARCTIS will be released under creative commons license, free of charge, via website (http://luftbildarchiv.univie.ac.at)

Using Low Resolution Satellite Imagery for Yield Prediction and Yield Anomaly Detection

Low resolution satellite imagery has been extensively used for crop monitoring and yield forecasting for over 30 years and plays an important role in a growing number of operational systems. The combination of their high temporal frequency with their extended geographical coverage generally associated with low costs per area unit makes these images a convenient choice at both national and regional scales. Several qualitative and quantitative approaches can be clearly distinguished, going from the use of low resolution satellite imagery as the main predictor of final crop yield to complex crop growth models where remote sensing-derived indicators play different roles, depending on the nature of the model and on the availability of data measured on the ground. Vegetation performance anomaly detection with low resolution images continues to be a fundamental component of early warning and drought monitoring systems at the regional scale. For applications at more detailed scales, the limitations created by the mixed nature of low resolution pixels are being progressively reduced by the higher resolution offered by new sensors, while the continuity of existing systems remains crucial for ensuring the availability of long time series as needed by the majority of the yield prediction methods used today.JRC.H.4-Monitoring Agricultural Resource

Estimation of Leaf Area Index Using DEIMOS-1 Data: Application and Transferability of a Semi-Empirical Relationship between two Agricultural Areas

This work evaluates different procedures for the application of a semi-empirical model to derive time-series of Leaf Area Index (LAI) maps in operation frameworks. For demonstration, multi-temporal observations of DEIMOS-1 satellite sensor data were used. The datasets were acquired during the 2012 growing season over two agricultural regions in Southern Italy and Eastern Austria (eight and five multi-temporal acquisitions, respectively). Contemporaneous field estimates of LAI (74 and 55 measurements, respectively) were collected using an indirect method (LAI-2000) over a range of LAI values and crop types. The atmospherically corrected reflectance in red and near-infrared spectral bands was used to calculate the Weighted Difference Vegetation Index (WDVI) and to establish a relationship between LAI and WDVI based on the CLAIR model. Bootstrapping approaches were used to validate the models and to calculate the Root Mean Square Error (RMSE) and the coefficient of determination (R2) between measured and predicted LAI, as well as corresponding confidence intervals. The most suitable approach, which at the same time had the minimum requirements for fieldwork, resulted in a RMSE of 0.407 and R2 of 0.88 for Italy and a RMSE of 0.86 and R2 of 0.64 for the Austrian test site. Considering this procedure, we also evaluated the transferability of the local CLAIR model parameters between the two test sites observing no significant decrease in estimation accuracies. Additionally, we investigated two other statistical methods to estimate LAI based on: (a) Support Vector Machine (SVM) and (b) Random Forest (RF) regressions. Though the accuracy was comparable to the CLAIR model for each test site, we observed severe limitations in the transferability of these statistical methods between test sites with an increase in RMSE up to 24.5% for RF and 38.9% for SVM

Estimating dry biomass and plant nitrogen concentration in pre-Alpine grasslands with low-cost UAS-borne multispectral data – a comparison of sensors, algorithms, and predictor sets

Grasslands are an important part of pre-Alpine and Alpine landscapes. Despite the economic value and the significant role of grasslands in carbon and nitrogen (N) cycling, spatially explicit information on grassland biomass and quality is rarely available. Remotely sensed data from unmanned aircraft systems (UASs) and satellites might be an option to overcome this gap. Our study aims to investigate the potential of low-cost UAS-based multispectral sensors for estimating above-ground biomass (dry matter, DM) and plant N concentration. In our analysis, we compared two different sensors (Parrot Sequoia, SEQ; MicaSense RedEdge-M, REM), three statistical models (linear model; random forests, RFs; gradient-boosting machines, GBMs), and six predictor sets (i.e. different combinations of raw reflectance, vegetation indices, and canopy height). Canopy height information can be derived from UAS sensors but was not available in our study. Therefore, we tested the added value of this structural information with in situ measured bulk canopy height data. A combined field sampling and flight campaign was conducted in April 2018 at different grassland sites in southern Germany to obtain in situ and the corresponding spectral data. The hyper-parameters of the two machine learning (ML) approaches (RF, GBM) were optimized, and all model setups were run with a 6-fold cross-validation. Linear models were characterized by very low statistical performance measures, thus were not suitable to estimate DM and plant N concentration using UAS data. The non-linear ML algorithms showed an acceptable regression performance for all sensor–predictor set combinations with average (avg; cross-validated, cv) R2cv of 0.48, RMSEcv,avg of 53.0 g m2, and rRMSEcv,avg (relative) of 15.9 % for DM and with R2cv, avg of 0.40, RMSEcv,avg of 0.48 wt %, and rRMSEcv, avg of 15.2 % for plant N concentration estimation. The optimal combination of sensors, ML algorithms, and predictor sets notably improved the model performance. The best model performance for the estimation of DM (R2cv=0.67, RMSEcv=41.9 g m2, rRMSEcv=12.6 %) was achieved with an RF model that utilizes all possible predictors and REM sensor data. The best model for plant N concentration was a combination of an RF model with all predictors and SEQ sensor data (R2cv=0.47, RMSEcv=0.45 wt %, rRMSEcv=14.2 %). DM models with the spectral input of REM performed significantly better than those with SEQ data, while for N concentration models, it was the other way round. The choice of predictors was most influential on model performance, while the effect of the chosen ML algorithm was generally lower. The addition of canopy height to the spectral data in the predictor set significantly improved the DM models. In our study, calibrating the ML algorithm improved the model performance substantially, which shows the importance of this step